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pJL^ ■ Abstract. Because the time needed for a simulation in lattice QCD varies at a rate exceeding 

the fourth power of the lattice size, it is important to understand how small one can make a 
lattice without altering the physics beyond recognition. It is common to use a rule of thumb 
that the pion mass times the lattice size should be greater than (ideally much greater than) four 
' (i.e., m^L S> 4). By considering a relatively simple chiral quark model we are led to suggest that 

£N| , a more realistic constraint would be m^(L — 27?) S> 4, where R is the radius of the confinement 

■ region, which for these purposes could be taken to be around 0.8-1.0 fm. Within the model we 

' demonstrate that violating the second condition can lead to unphysical behaviour of hadronic 

properties as a function of pion mass. In particular, the axial charge of the nucleon is found to 
decrease quite rapidly as the chiral limit is approached. 
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1. Introduction 



Our current capacity to compute hadron properties in lattice QCD is constrained by the need 
to take many limits. For example, we need to take the spacing a —* 0, the size of the lattice 
L — > oo and the quark mass to around 5 MeV. Of course, these limits are not unconnected 
because the size of a region of space big enough to contain (say) a proton will need to grow 
with the Compton wavelength of the pion. At present, with lattice spacings that represent 
a reasonable approximation to the continuum limit and for full QCD with reasonable chiral 
symmetry, the time for a lattice simulation scales roughly like m~ 9 and we are limited to pion 
masses larger than 0.4-0.5 GeV. This has led to much effort to explore the application of chiral 
perturbation theory as a tool for extrapolating hadron properties to the physical pion mass [1—3]. 

In such an environment there is great interest in seeing whether one can lower the pion mass 
without increasing the size of the lattice, thereby saving a factor of m\. Indeed, there have been 
many calculations for which even the rather optimistic rule of thumb that m n L > 4 has not 
been satisfied. Considerable attention has been devoted to studying such systems within the 
framework of effective field theory in order to understand how the relevant path integral might 
change as we go from one regime to another [4,5]. 

The question we ask is somewhat different. We consider a simple chiral quark model upon 
which we impose boundary conditions which roughly approximate those on a lattice. Simple 
inspection of the solutions naturally leads one to conclude that the condition noted above is 
incorrect. Indeed, the pion cloud of the nucleon does not even begin until one is outside 
the region of space in which the valence quarks are confined. Within the bag model this is 
characterised by the bag radius, R, and within the cloudy bag model as well as the model 



considered here, this radius is where the pion field peaks. The asymptotic behaviour of the pion 
field is therefore <j>(r) ~ exp[— m w (r — R)]/r 2 and the correct condition for the pion field to be 
small on a spherical boundary surrounding the nucleon is that (L/2 — R) S> ra^ 1 , with L the 
diameter of the spherical "lattice". (Note that in sects. 2 and 3 we will use L to denote the 
radius of such a region.) 

If we consider a typical case where the pion mass is large, say m w = 600 MeV, then on a 2 fm 
lattice m n L ~ 6 and the commonly quoted condition would suggest that we had a sufficiently 
large lattice. However, with a bag radius of order 0.8 fm, we find m n (L/2 — R) ~ 0.6 and the pion 
field has no chance to drop to zero before we reach the edge of the lattice! Indeed, even at this 
relatively large mass a 3 f m lattice would be a minimal requirement for full QCD simulations. 
From the mathematical point of view, having the boundary too close means that we are not 
restricted to the well behaved solution of the second order differential equation but can have a 
significant coefficient for the divergent solution. Through the coupling to the confined valence 
quarks this can in turn change the internal valence structure of the hadron. 

Considerations such as these help us to understand why some hadronic properties exhibit a 
dramatic volume dependence as the pion mass is varied. Perhaps the most famous example is 
the axial charge of the nucleon where early simulations revealed a striking decrease of qa as the 
pion mass decreased — in the opposite direction from the experimental data [6-8]. We shall see 
that our simple chiral model is able to reproduce this feature. 

2. A Simple Chiral Quark Model 

One of the earliest attempts to restore chiral symmetry to the MIT bag model was made in 
1975 by Chodos and Thorn [9]. Their solution, a precursor to the CBM [10,11] was a straight 
forward generalization of the linear sigma model, with pion and sigma fields coupling linearly 
to the confined quarks at the bag boundary. As well as presenting a perturbative solution in 
the pion field (as in the CBM), Chodos and Thorn attempted to find an exact solution to the 
resulting equations of motion. The only case where this was feasible was for a highly idealized 
baryon called the "hedgehog" , in which the 3 confined quarks all inhabit the same mixed spin- 
flavour singlet, and which hence is not an eigenstate of either isospin or angular momentum. As 
a result, the hedgehog solution does not correspond to any physical particle. Despite this fact, 
the hedgehog approximation does provide a simple classical solution with which we can then 
study a number of phenomena with relative ease. 

Chodos and Thorn suggest the following chirally invariant Lagrangian density based on the 
MIT bag model Lagrangian: 

C CT = $ijhl> - B]0 V - Xi;(a + if ■ 7f 75 )^5 5 + 1(^,7) (c^V) + 1(^7?) • (c^tt), (1) 

where ip, n and a are the quark, pion and sigma fields, respectively, and A is a Lagrange multiplier 
which turns out to be \(a 2 + ir 2 )~ 1 / 2 . This is invariant under the appropriate infinitesimal chiral 
transformation and the corresponding axial current: 

A* = -tpfj^ipOv + (0"<t)7t - a(d»Tr), (2) 

is conserved. 

To this Lagrangian we will also add explicit chiral symmetry breaking quark and pion 
mass terms. To generate a pion mass we could follow the method of the linear sigma model, 
spontaneously breaking chiral symmetry and then tipping the Mexican hat potential, generating 
masses for both the tt and a fields. However, this process leads to some very complex equations 
of motion for the pion and sigma fields, destroying the simplicity of Chodos and Thorn's classical 



solution. For simplicity we include the pion mass (and the corresponding quark mass) by hand, 
leaving the a field massless. The resulting Lagrangian is 



m q )ip — B]9v — Xip(a + ir ■ ttj^i/} 5s 
+\{d^Wa) + \{d^) ■ (^t?) - \m 2 J ■ vf. (3) 



Minimizing the action leads to five Euler-Lagrange equations. For a static spherical bag of 
radius R and static a and tt fields they are as follows: 

(t£ - m q )i/) = 0, r < R; (4) 
if ■ jtp = — £(a + ir ■ TTj^)ip, r = R; (5) 

V 2 ct = -^5(r-R); (6) 
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V 2 7? - m 2 7? = -£,ipifj 5 ^5(r - R); (7) 



1 d 

B = -^Q-ma + iT-7T 75 )^] r=Rt (8) 



where we use the notation £ = [cr 2 (i?) + 7r 2 (i?)] -1 / 2 . The first of these equations is just the 
Dirac equation for the confined quarks, whilst Eqs. © and © are equations of motion for the 
sigma and pion fields. There are also two boundary conditions — the linear boundary condition, 
Eq. and the non-linear boundary condition, Eq. JSJ), with similar roles to their roles in the 
MIT model. 

Chodos and Thorn solve this system of equations (without quark and pion masses) on an 
infinite volume. Following their method, we will solve the system on a finite spherical volume of 
radius L. Ideally we would use a rectangular volume similar to that used in lattice simulations 
but because of the radial nature of the solutions it is natural to use a spherical volume and any 
attempt to solve on a rectangular volume would be significantly more complicated. 

The quark field solution is just the MIT quark wavefunction, 



( ■ f nr 



Qr R \ Xh 0(R-r)e- iat / R , (9) 



where we define a = ER and 0, = kR (where k is to be determined by the appropriate eigenvalue 
condition) related by the equation 



a = ^Q 2 + (m q R) 2 , (10) 

and 

a ± m n R 

a ± = J q —. 11 

V a 

The defining property of the hedgehog solution is the choice of static, radially dependent 
sigma and pion fields, 

n(r) = g(r)f, (12) 
a(r) = f(r), (13) 

and the choice of spinor-isospinor, \h, proportional to a mixed spin-isospin singlet state, 

x k = -^(|dT>-M», (14) 



so that it has the property, 

(<? + f) Xh = 0. (15) 

These choices greatly simplify the equation for the pion field by ensuring that ^rj^ip is 
proportional to f and hence to tt: 

TpTj 5 iJj = -2ij ji xlxh r. (16) 

Substituting Eqs. @, (JT2|) and (fT3)) into the equations of motion for 7? and a, yields the following 
two second order differential equations in / and g: 

f"(r) + -f'(r) = b(alj 2 (n)-a 2 j 2 l (n))xlx l J(R-r) 

= a5{R-r), (17) 

g"(r) + -g'{r)-(\ + ml)g{r) = £a + a_j (n)ji(fi)xtx/.*(fl- 

= bd(R-r). (18) 

Enforcing regularity at the origin, Eqs. (|17|) and (|18|) have solutions 

/(r) = f + a R 2 (^--)e(r-R), (19) 
R r 



cosh rm, sinh rm„ 



g (r) = ( - + — ^r^) C + b(R sinh i?m, -)0(r - R) 

rm^ r z m z 



cosh i?m„ 

m„ 



+ 



( Sinhrm - - C ° S 9 h 7-)b(^cosh^ - Sinh ^)g(r - R), (20) 
rm, Hm/ m, 

where /o and C are constants of integration. 

At the boundary r = L we set the derivative tt' to zero to resemble the periodic boundary 
conditions applied in lattice simulations. We stress that in this context L is not the lattice 
length, which would be of order 2L. Because there is no sigma mass we cannot enforce a similar 
condition on a. Instead we leave /o as a free parameter for the moment. The resulting solutions 
are 

f(r) = f o + a R 2 (^--)0(r-R) 1 (21) 
R r 

sinh Rm n sinh rm, coshrm 

g{r) = b(RcoshRm n )( ^ — ^ — )8(r — R) + 

m, rm^ r z m z 

un ■ l. n cosh i?m cosh rm, sinhrm 

b(RsmhRm n )( ^ — ^ — )"{R — r) + 

m n rm n r z m z 

sinh flm cosh rm, sinhrm, 

b(RcoshRm^ )( 1 5 — i Lm, , (22) 

rm, r z m, 

where the L-dependence is given by the function, 



cosh Lm w n sinh Lm n • o cosh Lm, 
z — r2™2 r 



r(LmJ = -^f Vr 2 • (23) 

I sinh Lra, o cosh Lra, , osinhLm, / 
— z — r2™2 r ^ 



Lm, L 2 m^ LPm% 

nh Lm, 2 cosh Lm, 1 2 sinh Ln 

(Note that r(Im.) — ► 1 as L — > oo). 



2.1. The linear boundary condition 

The solutions to Eqs. ((3J), © and are dependent on the parameters f2, R, L, m q and but 
also on the unknown values f(R) and g(R) = go, which are given by 



f(R) = fa, 

g{R) = -£,Oi + a_j Q jix{xh{RcoshRm 1T - 
sinh Rm n cosh Rm^ 

= -£d(Q,m q ,m n ,R,L). 



sinh /2m.,, 

cosh i?m^ 
Rm^ 



x 

sinh i?m 



R 2 m? 



-)r(Lm„ 



(24) 

(25) 
(26) 



If we now apply the linear boundary condition, Eq. JSJ), we may express /q , g 2 . and l/^ 2 in terms 
of y and d: 



fa 



±y l 



d, 



i-y z 
i + y 2 



d. 



l 



i + y z 
i-y 2 



d, 



(27) 



so that our solution now depends only on the parameters f2, m„, m^, R and L. 



2.2. The non-linear boundary condition 

The last condition to apply is the non- linear boundary condition, Eq. @, which sets the bag 
energy density B. Because of the discontinuity in the derivatives of the tt and a fields at the 
bag boundary, r = R, we take the average of the derivative on each side 1 , that is: 



da 
dr 



1 /da^ do-(° ut 1 
- 1 — - — + 



r=R 



r=R 2\dr dr 
With the solutions for tp, jf and a, Eq. Q leads to an expression for B of the form, 

A-kR a B = t(Q, m q , m^, R, L), 



(28) 



(29) 



where r is a complicated function of the bag frequency f2, hedgehog and volume radii {R and L 
respectively), and quark and pion masses. Because of its complexity, the explicit form of r and 
the involved steps required to find it are not shown. 

Once the masses and volume size are set, Eq. (|29p becomes a three-way relationship between 
the parameters B, R and SI. By choosing a value for the bag energy density (a property of 
the vacuum) the radius R, and hence the entire solution, will be completely determined by f2. 
Using this fact we can now find an eigenvalue condition for O. Using the linear and non-linear 
boundary conditions, we find that the value of the sigma field at the boundary r = L, 



1 1 



a{L)=f(L)=h + aR\-- 1 ), 



(30) 



is completely specified by the energy parameter Q. Fig. ^shows a plot of 1//(L) 2 against Q 
for fixed parameters B, m. 



m n and L. Physical solutions occur for Q within specific allowed 
regions or bands. Outside of these regions either a 2 is negative and hence a is imaginary, or 
t(S7) is negative, in which case there are no real solutions to the n.l.b.c. By selecting a boundary 
condition for the sigma field (represented by a horizontal line in Fig. we can find eigenvalues 
of the bag frequency, 0. From Fig. [I]we see that l/f(L) 2 has a limited size so there is a limit 
to how small we can make o~(L). 



1 With this prescription, Chodos and Thorn found that the non-linear boundary condition (n.l.b.c. ) corresponded 
with conservation of energy in the hedgehog. 
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Figure 1. Plot of the inverse square of the value of the sigma field at the boundary r = L 
versus the bag frequency f2. The horizontal line represents the boundary condition f(L) = f^. 



So what boundary condition should we apply to cr? Starting with the PC AC relation 



we substitute the axial current of Eq. (J2J into the left-hand side and take the limit r — > oo 
(where ip = 0) then we find 



Comparing with Eq. (|31[) . we see that the vacuum expectation value of the sigma field must be 
(a) = f n so we expect the sigma field to go to at infinity. 

To find the hedgehog solution on an infinite volume we would therefore choose /o such that 
/( r ) ~~ * /x as r — > oo. For large volumes, setting f(L) = f n should still be an appropriate 
boundary condition. But for smaller volumes this is too harsh a condition and is found to 
introduce large volume dependence. As an alternative, we fix the sigma field at r = L to the 
value taken by cr(L) in the infinite volume solution (in which a — > f„ as r — > oo). This results 
in a more modest volume dependence, which seems appropriate given that in non-linear chiral 
models the sigma mass is usually set to infinity and we expect volume dependence to arise 
mainly from the effect of the tt field. 

The solution displayed in Fig. [2] for the quark ground state and associated tt and a fields was 
found using the latter boundary condition on a volume of radius L = 100 fm. The quark and pion 
masses were set to m q = 5 MeV and m w = 140 MeV and the energy density B = 13.27 MeV/fm 3 
yields a bag radius of R = 1.4 fm as L goes to infinity. In this case the lowest energy state is 
given by the eigenvalue Vt = 1.588. Although the choice of bag radius R = 1.4 fm may appear 
slightly large, this is in fact the smallest value for which ground state solutions can be found on 
an infinite volume in this model. Clearly what matters in terms of physics consequences is the 
difference between L and R, rather than the absolute values. 

2. 3. The mass of the hedgehog 

The energy of the combined system is found by integrating the energy-momentum tensor and 
can be divided into three separate pieces, 



(0|VMk(<z)} = Lm 2 J ab e 



(31) 



(oivTOM?)) 



(0\{(d 2 a)7r a -a(d 2 7r a )}\7r b (q)) 



(32) 
(33) 



Ehh — + Sbag + -Ew- 



(34) 




Figure 2. The quark density i/j^i/j (full line) and magnitudes of the pion (dashed line) and sigma 
(dotted line) fields in the hedgehog solution calculated at the physical pion mass on a spherical 
volume of radius L = 100 fm. 
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The quark and bag energies we already know from the MIT model: 

(35) 
(36) 

The pion and sigma energy is given by the integral 

E„ = \f d\[(Va) 2 + (VK) 2 + m 2 J 2 ] (37) 

= 2vr f dr[r 2 f'{rf +r 2 g'{rf + {2 + r 2 ml)g(rf\. (38) 
Jo 

For the solution presented in Fig. |2 the total energy (or hedgehog mass) is 998 MeV, which is 
comparable to the physical nucleon mass M N ~ 940 MeV. 

The axial coupling constant is found by integrating this nucleon matrix element of the axial 
current over all space. In the hedgehog we integrate A^ h {x) (given by Eq. @): 

g\ h 5ij = j d\(hh\A){r)\hh) (39) 

tfV[^-757V + i&<r)*3 ~ ( 4 °) 

Applying the hedgehog solutions, one can obtain closed forms for the quark and meson 
contributions. 



3. Results: Volume Dependence of the Hedgehog 

Using the technique described in the previous section, we can find solutions for the groundstate 
hedgehog on a range of different sized spherical volumes and at different pion and quark masses. 
For each solution we calculate the hedgehog energy and the axial coupling constant. The bag 
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Figure 3. The quark density ^i/j (full line) and magnitudes of the pion (dashed line) and sigma 
(dotted line) fields in the hedgehog solution calculated at the physical pion mass on spherical 
volumes of radius 2.5 and 3fm. 



energy density is fixed at B = 13.27 MeVfm" 3 (to give R = 1.4 fm in the infinite volume solution) 
and the quark and pion masses fixed in the proportionality defined by the Gell-Mann-Oakes- 
Renner relation so that m q = (m^-/140MeV) 2 5MeV. By choosing the volume size, L, and pion 
mass, m„, and applying the boundary condition on a we find the lowest eigenvalue of f2, which 
then completely specifies the solution. 

At small pion masses, if the volume becomes too small then there is no groundstate solution 
which satisfies the boundary condition. As L decreases, the sigma field will grow whilst the 
value of <j(L) required by the boundary condition gets smaller, such that eventually there is no 
value of f2 for which cr(L) is small enough. At the physical pion mass we cannot find solutions 
on volumes much smaller than L = 2.5 fm. For larger pion masses the pion and sigma fields are 
smaller and we do not encounter this problem. 

Graphs of the solutions for 2.5 and 3fm at the physical pion mass are shown in Fig. El 
We observe that the hedgehog solutions on different sized volumes are relatively similar. 
Interestingly, as the volume size decreases there is a small increase in the bag size, accompanied 
by a decrease in the quark density over the bag. There is very little change in the n and a fields 
inside the bag volume. The most noticeable change is in the behaviour of the pion field outside 
the baryon, which becomes much larger for smaller volumes because there is less distance over 



which it can flatten out to satisfy the periodic boundary condition. 

The results for the axial coupling constant g A are plotted against m\ in Fig. 0J with each line 
representing a different volume. We see that the results for different volumes converge at large 




Figure 4. Plot of the axial coupling constant g A in the hedgehog versus ml, as calculated 
on spherical volumes of radius 10 fm (open star), 5fm (filled star), 3fm (open triangle), 2.5 fm 
(filled triangle), 2fm (open diamond) and 1.8 fm (filled diamond). 

ml, revealing very limited volume dependence in this region. This is what we expect to see 
because at large pion masses the pion field will die out quickly so that the effect of the boundary 
is minimal. At low the axial coupling constant exhibits large finite volume effects. Indeed, 
for the two lightest masses, = 140 MeV and m^ = 200 MeV, the effect is very large, with g A 
decreasing by around 30% and 50% from L = 10 to 2.5 fm in each case. It is also interesting 
to note the turn-around in the low m n behaviour of g A between the L = 10 fm and L = 5fm 
solutions. 

4. Discussion 

Our plot of the hedgehog axial coupling constant reveals a large volume dependence in the low 
m^ region, with g A at the physical pion mass decreasing by over 50% as the volume size decreases. 
This behaviour is similar to that observed in recent lattice QCD results. From the plots of the 
hedgehog solution in Fig. El we see that the important region to consider is the distance L/2 — R 
(where we once again use L to denote the length of the side of the lattice) . This is the distance 
between the edge of the baryon and the volume boundary. Inside the baryon bag the fields do 
not vary much with volume size. We expect that this conclusion is far more general than the 
particular model considered. 

Clearly this has been a very simple study, involving a number of major approximations which 
must be remembered when considering our results. But by employing such a simple model we 
have been able to generate results over a wide range of volumes and pion masses with relative 
ease. Because of the exact nature of the solutions, at this point it would be very easy to 
examine the volume and mass dependence of other hedgehog properties. It would clearly be 



very interesting to employ a spontaneous symmetry breaking mechanism to include the pion 
and sigma masses. 
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